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Abstract 

In the rest frame of a many-body system, used in the calculation of its static and scattering 
properties, the center of mass of a two-body subsystem is allowed to drift. We show, in a model 
independent way, that drift corrections to the nucleon-nucleon potential are relatively large and 
arise from both one- and two-pion exchange processes. As far as chiral symmetry is concerned, 
corrections to these processes begin respectively at 0(q 2 ) and 0(q 4 ). The two-pion exchange 
interaction also yields a new spin structure, that promotes the presence of P waves in trinuclei 
and is associated with profile functions which do not coincide with neither central nor spin-orbit 
ones. In principle, the new spin terms should be smaller than the 0{q^) spin-orbit components. 
However, in the isospin even channel, a large contribution reverts this expectation and gives rise 
to the prediction of important drift effects. 
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I. INTRODUCTION 



This work is motivated by a private question posed by Alejandro Kievsky some years 
ago, concerning the possibility of novel forms of spin dependence in the nucleon-nucleon 
interaction, when one is not in the center of mass frame of the two-body system. In the 
study of static and scattering properties of many-body nuclei, calculations are performed 
in the rest frame of the larger system and a two-body subsystem is allowed to drift. This 
picture led him to introduce a phenomenological three body forcefl], which improved the 
description of the N — d vector analyzing power A y . 

Nowadays, the outer layers of the NN interaction, represented by one-pion and two-pion 
exchange potentials [OPEP and TPEP), are set in solid foundations due to the use of 
chiral symmetry. Nuclear processes are dominated by the light quarks u and d, and one is 
not far from the massless limit, in which QCD becomes invariant under both isospin and 
chiral SU(2) x SU(2) transformations. Chiral symmetry is realized in the Nambu-Goldstone 
mode and the QCD vacuum can bear collective excitations, identified as pions. A suitably 
formulated chiral perturbation theory (ChPT) allows deviations from the massless limit to 
be treated systematically [2]. As low-energy QCD calculations are prevented by its non- 
Abelian character, in practice one works with chiral effective theories, in which elementary 
nucleons interact by exchanging pions. 

In chiral perturbation one uses a typical scale q << 1 GeV, set by either pion four- 
momenta or nucleon three-momenta. The leading term[3] in the NN interaction is the 
OPEP, at O(q ). The TPEP begins at 0(q 2 ) and two independent expansions up to 
0(q A ) are presently available. One of them, based on heavy baryon ChPT [4], uses non- 
relativistic lagrangians from the very beginning and the inverse of the nucleon mass as 
an expansion parameter. The other one, produced recently by our group[5, 6], is based 
on relativistic expressions, written in terms of observable coefficients and covariant loop 
integrals. The use of a relativistic language frees one from particular reference frames and 
allows a straightforward treatment of two-body interactions in which the center of mass is 
able to move. Here, we rely on our previous work in order to derive the drift contributions 
to the NN potential. For the sake of definiteness, we stay in the realm of three-body nuclei, 
but results can be easily generalized to larger systems. 

Our presentation is organized as follows. In section II, we review the dynamical role of 
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two-body interactions in trinuclei. This sets the stage for the derivation of drift interactions, 
which is performed in section III. Results are summarized in section IV, whereas technical 
issues, concerning kinematics and spin operators, are left to appendices. 



II. DYNAMICS 

The interactions of a three-nucleon system in momentum space are represented by the 
operator W, defined by [7, 8] 

(p'i,P2,P' 3 \W\p 1 ,p 2 ,p 3 ) = - (2tt) 3 5 3 (p' 1 +p' 2 +p' 3 -p 1 -p 2 -p 3 ) t 3 (p' l ,p' 2 ,p' 3 ,p 1 ,p 2 ,p 3 ) , (1) 

where £3 is the proper part of the non-relativistic three-body transition matrix. In configu- 
ration space, the position of nucleon i is described by 7*j and one uses the Jacobi variables 



R= (ri + r 2 + r 3 )/3 
which correspond to 



r = r 2 -r 1 



p = (2r 3 -r 1 -r 2 )/V3. 



(2) 



P P P P p p 



P3 = y + 7f 



(3) 



The Schrodinger equation for the internal degrees of freedom is obtained by using P 
P' = and given by 



with 



m m 1 



V3' 



J drdpW(r',p';r,p) ip{r,p) , (4) 



W(r',p';r,p) 



(2n) 



12 



dQ r dQ p dq r dq p 



xe 



[Q r <r'-r)+ Q p -(p>-p) + q r -(r> + r)/2 + q p (p> + p)/2] ^ Q ^ q 



pi 1 



(5) 



and Qi = and q { = (p--pj, for i = (r,p). 

In this work we are interested in describing two-body interactions between nucleons 1 
and 2 and note that the conservation of p 3 implies q p = 0. We write 



t3(p'l,P2,P3iPliP2iP 3 ) = W 



s 3 (q P )HQ P ,QnQr) , 



(6) 



where t 2 is the two-body t-matrix, and the corresponding potential becomes 



W 2 (r',p';r,p 



dQ r dQ p dq r 



xe 



(2vr) 9 LV3. 
Q r -(r'-r» Q p -(p>- P ) + q A r> + r)/2} Q g j 



(7) 



In isospin space, the amplitude t 2 reads 



(8) 



The usual spin decomposition is obtained by going to the center of mass frame of the 
two-body system, where one finds 



4] 



LS ,± 



rrr 



ss ,± 



(9) 



with two-component operators defined by 



n LS = i{(T^ + ( T^)-q r xQ r l2 



n T = - q l (3 *^.q r *M.q r - 
n Q = 4aW-q r xQ r <r®-q r xQ r . 



(10) 

(11) 
(12) 

(13) 



In this formulation, the two-body interaction does not depend on Q p and is completely 
decoupled from the larger system it is immersed in. The Fourier transform of this result 
produces the configuration space potential, given by 



W 2 (r\p';r,p)=S 6 (r'-r) 5 A (p'-p) 



V3. 



cm = V c + Vts ^ls + V* Qss + V± Q T 



(14) 
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where we have kept only local and spin-orbit contributions and the spin operators read 



n LS = L-(<T^+*M)/2, 

^S = <X (1) -<T (2) , 



The radial functions are given by 



(15) 
(16) 
(17) 



with x = fir and 



Uf{x) 



V±(r) = U±{x) 



V±(r) 



d 3 q 



!*l±u± ( X ) 

m 2 x dx LSK ' 



2 d 

oc doc 



dx 2 

d 2 1 d 

dx 2 x dx 



UssW 



U$(x) 



(18) 
(19) 

(20) 
(21) 



I = {C,LS, SS,T}. 



(22) 



As we discuss in the sequence, the fact that the two-body CM is allowed to drift gives 
rise to extra interaction operators in the potential. 



III. DRIFT TERMS 



Corrections to the NN potential due to the motion of the CM are derived by evaluating 
T, the covariant t-matrix for the on-shell process N(pi) N(jp2) ^(p'i) N{p' 2 )-i and writing 
the result in terms of two-component spinors, using the expressions of appendix B. This 
gives rise to an amplitude expanded in terms of Pauli spin operators. Dividing it by the 
factor 4mE present in the relativistic normalization, one obtains the amplitude i 2 , which is 
to be fed into eq.(7). In this work we concentrate on contributions from processes due to 
the exchanges of one and two pions. 

The transformation of a t-matrix into a potential to be used in a dynamical equation is 
not trivial and depends on a number of important conventions. These range from the very 
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nature of of the equation adopted to tacit assumptions concerning the off-shell behavior of 
the potential. The latter class of effects appears as corrections to leading order effects and 
was discussed in a comprehensive paper by Friar [9]. We here stick to the conventions used 
long ago by Partovi and Lomon[10, 11]. 

• OPEP 

The covariant amplitude for on-shell nucleons reads 

T = rW- [ii75ti ]<i) [u 75 up) , (23) 

ft T ~ V 

where g^, f n , fi>, m are respectively the axial and pion decay constants, the pion and nucleon 
masses. Using eq.(A8) for the momentum q and eq.(B4) for the spinor matrix element, one 
finds the two-component amplitude 

r = (l). (2) 9 2 A ™ 2 1 

ft {ql+^)-A{q r -Q p y/ZE* 

M 2 [[m+S/2 + 2Q r -Q p /V3S] a^-q r - 2 q r Q p tr™ ■ (Q r + Q p /V3)/V3 s} 
{[m+S/2-2Q r -Q p /V3S} ( r^-q r + 2q r -Q p a^-(Q r -Q p /V3)/V3S} , (24) 
where M 2 is a normalization factor, 

N 2 = {[(m+£/2) 2 + [A(Q r -Q p ) 2 -(q r -Q p ) 2 ]/3S 2 ] 2 - 16(m+£ /2) 2 (Q r ,Q p ) 2 /SS 2 }^ 

(25) 

and £ is the total energy of the two-nucleon system, determined by the condition 

£ 4_ 4(m 2 + q2 + 4Q 2 + Q 2 /3) ^ +(4/3) [ {qr .Q p f + A{Qr .Q p f] = . (26) 

This t-matrix is fully relativistic and contains no approximations. All its terms involving 
the variable Q p vanish in the rest frame of the two-body system and therefore can be 
interpreted as drift effects. With the provisos discussed in ref. [9], it could already be used 
as input into a momentum space dynamical equation. Alternatively, in the framework of 
chiral perturbation theory, one might wish to rewrite it as a power series, truncated at a 
given order. 
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In configuration space, the variables Q r and Q p correspond to non-local operators and 
are usually associated with gradients acting on the wave function. In order to restrict the 
corresponding complications to a minimum, we remain in the limited scope of eq.(4) and 
keep only terms linear in these momenta. This procedure is referred to as the linear gradient 
approximation. 

Within this approximation, the OPEP retains its usual local form, given by 

_ _ Iga Q aTTI 1 t f c\l~r\ 

t ss = —t T ~ — — + local corrections . (27) 

• TPEP 

Quite generally, for each isospin channel, the spin content of the TPEP is given by [5] 



T* = [««]« [uup {T± D ) - ^ [««]« [ua, x {p'-pYu]V (T± B Y 

ua» x (p'-pyu]V [ua vp (p'-pyu]M (Zf B ) A " , (28) 



4m 2 



where the functions 1 involve loop integrals and have a Lorentz structure realized in terms 
of the kinematical variables W, z and q, defined in appendix A. Terms proportional to q do 
not contribute for on-shell nucleons and we have 



1 db) ~ 2m 2m ' [ZJ) 

1 bd) ~ 2m 2m ' ^ 

(T ± \*P - „V t(9)± + WXWP tH± + — tW± fS1 1 

V-BBi ~ 9 BB "I" _B_B T ^ ^BB ■ i 61 ) 

The amplitudes X were explicitly calculated in ref.[5], as functions of the invariants W 2 , 
z 2 and q 2 , and the two-pion exchange interaction is described by 
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[uu] m [uu]W 



<r± , J_ T H± i q 
- L DD ' " i nn ~r 



2m 



2 DB 



16m 4 



(w)± 



BB 



-—{[uu\V[u a, x (p'-pY u] W - [u a, x (p 1 ' -pf u\^[uu] (2) } 



2m 



T (w)± T (z)± _q_ T (w)± 



-DB 



4m 2 



4m 2 



^ [ua, x {p'-pYu\^ [ug vp {p'-pYu]W 



2 2 ^ + 4m2 I 



_7-0)± , <7-(*)±> 
X BB 



-BB 



(32) 



This result can be recast in a form similar to eq.(24), by using the spinor matrix elements 
given in appendix B. One finds 



T ± = N 2 [l± D [2m(m+S/2)-2(q r -Q p ) 2 /3£ 2 + q 2 /2 - i ct^ -q r x (Q r + Q p /V3) 
2m(m+S/2)-2(q r ,Q p ) 2 /3S 2 + q 2 ./2 - ict^ -qrXiQr-Qp/Vz)] + ■ ■ ■} , (33) 



and its full drift content becomes explicit. However, for the sake of simplicity, we remain 
in the framework of the linear gradient approximation. Using eqs.(A6-A8), one learns that, 
in this case, the variables W 2 , z 2 and q 2 do not depart from their CM values and the only 
sources of drift corrections are the spin functions. The results of appendix B yield the 
following non-relativistic amplitude 



tf 



m 



2 D ' 



where the drift operator fi^ is given by 



(34) 



(1)_«.(2)> 



■q r xQJ2V3 



(35) 



and the profile functions read 



D 



m \ 4m / qi 

7 — 1+ # 



T ± Qr T M± 

2m 2 



Qr T {w)± 
BB 



16m 4 



+ 



e Qr T (g)± 



m A 



2 ^BB 



(36) 



with A 2 = 4m(e+m) and e = \/ m 2 + qr 2 +4Q 2 .. This result is fully model independent, since 
it springs directly from Lorentz covariance and is constrained just by the linear gradient 
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approximation. The profile functions t% do not coincide with any other components of the 
TPEP which, in the same approximation, are given by[5] 



<tr T (w)± 
DB 




(w)± 



lQm^ BB 



Qr 



(s)± 



1 J_ 2(lr \ 



T 



16m^ BB 



(*)± 

DB 



(37) 



Am 2 



\ T {w)± 

' y2 I BB 



Am 2 



1+ 



Am? 



X 2 



1 



(s)± 

BB 



m 
e 



12 



■X 



0?)± 

BB 



(38) 



(39) 



We consider here the expansion of the TPEP to 0(q 4 ), using eqs.(9) and (34), 
which requires t@ — > 0(g 4 ) and t$,tg S , t^} — > C(g 2 ). The expansion of the 

various profile functions is performed using the results {X^ + ,X^ + ,X^ + } — > O(g ), 

{^ + ,^,^2" ^i" ^iTHOC? 1 ), {^}-0(s 2 ), &HO(? 3 ), {41",4t} ~ o, 

and one finds 



J D 



J D 



III 

e 

m 

e 



Vr T (9)+ , I 

im 2 + 2 



<lm 2±DB 



{0(q 2 )+[0(q 3 )}} 



-1 



DD 



{0(q 2 )} • 



(40) 



(41) 



In the expression for tj, the term within square brackets is 0(q 3 ). Nevertheless, we have 
kept it, for it is anomalously large. Considering comparable terms in eq.(37), one writes 



III 

e 



T + Qr T {w)+ 

2m 2 



171 ) 4_ T< x )+ 

e 



'DB 



9r ( T (w)+ 

Am 2 I ^ 2 ^ BB 



tc = -{l DD }^{0( q 2 )} 



t 



m 
e 

m 



LS 



1 



(«>)- 



DB 



+ T {z) ~ - 

' DB 



1 



{0(q 3 )} , (42) 

" Xn }-{0(g) + [0(g 2 )]} , (43) 

(44) 

{O(q) + [0(q 2 )]} . (45) 



DD 



These results show that the drift potential has little affinity with the spin orbit term and, 
at the chiral order considered, can be written as 
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3q2 i 

t D = J^J t ss - 2 l c ' ( 46 ) 



t D = ~\tc- (47) 
The Fourier transform of eq.(34) yields the configuration space structure 

V(r) ± =V(r) ± ] m + V$Sl D , (48) 



with 



^D = ^(<x (1) -<r (2) )-rxHV p ), (49) 

V^(r) = ^- 2 -^U^(x), (50) 

m z x ax 

U$<?) = -j 0js eiQr ' rt »(<lr). (51) 

The presence of the operator 

V p = V p - V p (52) 

in eq.(49) ensures that results are symmetric under the exchange of initial and final states. 
Using results (46) and (47), one has 



U£ = \v+ s - l -V+ , (53) 

U D = -\V C -. (54) 

In figs.l and 2 we display the profile functions for the drift and spin-orbit potentials, 
derived from our 0(q A ) expansion of the TPEP[5, 6]. These results do not include short 
range effects and cannot be trusted for r < lfm. We recall that both components of the 
force are multiplied by 0(q 2 ) spin operators and hence we need to keep just 0(q 2 ) terms 
in Vd- As shown in eqs. (40-45), in principle one should have V£ s ~ 0(q) > ~ 0{q 2 ). 
These expectations are confirmed in the figures, provided one uses the 0(q 2 ) dotted curve 
for . However, when the 0(q 3 ) term associated with the central potential is kept, one 
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has a complete subversion of the expected chiral hierarchy, associated with the prediction 
of a rather large effect in the isospin even channel. 

i 1 1 1 1 -i— — 1 — i 1 ==»». 1 0.0 



-10 - 



> 

CD 



-20 - 




- -0.1 



- -0.2 



12 3 r(fm) 4 

FIG. 1: (Color online) Isospin even drift (full and dotted lines) and spin-orbit (dashed line) poten- 
tials; the dotted line is 0(q 2 ) whereas the full one is 0(q 2 ) + 0(q 3 ). 



> 
CD 




0.00 



- -0.01 



-0.02 



FIG. 2: (Color online) Isospin odd drift (full line) and spin-orbit (dashed line) potentials. 
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In order to produce a feeling for the role of drift interactions in trinuclei, we note that 
their ground states contain S, P and D waves, but they are heavily dominated by the 
principal S component, which is fully symmetric under the exchange of nucleon coordinates. 
Using the notation of ref.[12], we write 

\S) = S(r,p)r 1 $, (55) 

where 

r^f = [\ m ~^)s - \m + fi)s |m~z')_r] (56) 

is the totally antisymmetric spin-isospin = (1/2, 1/2) wave function with third components 
fi and i, whereas \m + ) and \m~) represent respectively even and odd mixed symmetry states 
under permutation of particles 1 and 2. The leading term of the function S(r, p) is known[13] 
to depend just on the hyper-radius £ = \Jr 2 + p 2 and hence the most important coupling 
introduced by the drift potential is associated with the structure 



n D \s) ~ n D s® r$j* = ^ ^ {{- ® Y 1 (p)] 1 ® |m + >^ /2 

+V2 [[Ytir) (8) FxO)]! ® |s)s]£ /2 } «)/ + [[^(r) ® ^(p)], ® |m-) 5 ] A ; /2 |m" i)/}(57) 

1 3)5 being the spin 3/2 state. This result indicates that the drift potential enhances the role 
of P waves in trinuclei, as one might have guessed directly from eq.(49). 



IV. SUMMARY 



In nuclei containing three or more nucleons, the center of mass of a two-body subsystem 
is allowed to drift. This kind of movement does affect the forms of both one- and two-pion 
exchange contributions and gives rise to important non-local corrections to the potential. 
As interactions of this type are difficult to be dealt with in configuration space, we have 
restricted ourselves to the simplest possible non-local operators, proportional to single gra- 
dients acting on the wave function, which arise in two-pion processes. Using a relativistic 
chiral expansion of the two-pion exchange NN potential to 0(q 4 ) derived previously, we 
have shown, in a model independent way, that the profile functions of the drift corrections 
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do not coincide with none of its components. The spin dependence of the drift term is 
implemented by the operator 



1 



( -(l)_ cr (2)). rx (_ zVp ) , 



4v/3 



where r and p are Jacobi coordinates associated with two and three bodies. This structure 
promotes couplings between S and P waves, enhancing the role of the latter in trinuclei. 

As far as chiral symmetry is concerned, drift corrections begin at 0(q A ) and, in principle, 
should be smaller than spin-orbit terms, which begin at 0(q 3 ). However, in the isospin even 
channel, the same dynamical contribution that makes the its 0(q 3 ) central component to 
be larger than the 0(q 2 ) odd counterpart subverts the expected chiral hierarchy and gives 
rise to the prediction of important drift effects. 
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APPENDIX A: KINEMATICS 

The conventions used here are the same as in ref.[5]. The initial and final nucleon mo- 
menta are denoted by p and p' and we define the variables 



The interacting nucleons are assumed to be on shell and the following constraints hold 



W = Pl + P2 = P[ + P' 2 , 
Z= [(Pl+P'l)-(P2+P' 2 )}/ 2 

q = p[-pi=P2- p' 2 , 



(Al) 
(A2) 



(A3) 



m 2 = (W 2 + z 2 + q 2 )/A , 



(A4) 



W-z = W-q — z-q — . 



(A5) 
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Using the Jacobi variables defined in eq.(3), one has 



W = (S, -2Q p / v / 3) , (A6) 

z = (AQ r -Q p /SV3,-2Q r ) , (A7) 

q=(2q r -Q p /SV3,-q r ) . (A8) 

where £ is the total energy of the two-body system. If there were no drift, this energy would 
be written in terms of the single particle CM energy e as 

£ cm = 2e = 2^/m 2 + q 2 r +4Q 2 r . (A9) 

Explicit calculation yields 

£ 2 + (4/3) [{q r -Q p ) 2 /£ 2 + ^Qr-Q P Y/£ 2 -Q^ =4e 2 (A10) 
and hence, in the linear gradient approximation, 

£~2e. (All) 



APPENDIX B: SPIN OPERATORS 

We present here the changes induced in the spin operators due the drift of the two-body 
CM. With our conventions, we write 



[uTup = jA/V [E'+m,-a-p'] 
M = l/y/(E'+m)(E+m) 









E+m 




ap'] 






CT-p 


«) 



(Bl) 

;b2) 



for a generic Dirac matrix V. We display results for nucleon 1 and those corresponding to 
nucleon 2 are obtained by making q r — > —q r and Q r — > Q r . For the normalization, one has 



+£/2) 2 + A(m+£/2)Q r -Q p /V3S+[A(Q r -Q p ) 2 -(q r -Q p ) 2 }/3£' 



-1/2 



Iga , , v 

~ l/(m+e) , 



(B3) 
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where the last equality corresponds to the linear gradient approximation. 
The OPEP, eq.(23), is based on the function 



[u l5 u}^ = Af X ^{[m+£/2 + 2Q r -Q p /V3S] aW-q r -2q r -Q p tr<»-(Q r + Q fi /V3)/V3£} 



l M3 u-7 - ^tr^(i) 



X 
B4) 



The expression for the TPEP is given by eq.(31) and employs the operators 



p') ti(p)]W = {a/V [2m(m+S/2) - 2(q r -Q p ) 2 /3S 2 + q 2 r /2 - z a-q r x(Q r + Q p /V3)] x}^ 



u(p) u(p)} 



(1) _ J y t 
cm l A 



fe+m) 



■q r xQ /V3 



(i) 



(B5) 



[— u{p') a, (p'- P ru(p)}V = {{M/2m)x ] [(m+£/2)[q 2 - 2i aq T x (Q r+ Q p /V3~)} 



-2{q r -Q P ) qAQ r + QjV3)/V3£ 



X 



X' 



(1) 



— crq r xQ /Vz 



(i) 



(B6) 



[^-u(p') a M (p'-pY <P)} (1) = {^X ] Um+E/2) icrxq r + [-q 2 r + 2i a q r x (Q r + Q p /VS) 



2m 



+4(q r -Q p ) 2 /3£ 2 ](Q r + Q p /V3)/2m-(q r -Q p ) [q r + 2icrx(Q r +Q p /V3)}/V3£ 



X 



(i) 



q 2 r Q p /V3m + (q r -Q p )q r /V3e X \ ( B7 ) 



(i) 



These results allow one to write 



{[t2«P>[t2 u]< 3 >} -{■•■} 



4m 



(e+m) (e+m)' 



D 



■^[fi ^o/ - P r ^ ^p'-pt «] (2) ) ^ - {■ • -} cm 9 Xp - m ^ Tm) 



4m 2 



[u VAP-PY a up (p'-pY H (2) 



Z 7lP Iga 



Am 2 



Am 2 



(B8) 
(B9) 
WO) 
(Bll) 
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where the functions {• • -} cm are given by eqs.(A32-A35) of ref.[5] and the two-component 
spin operators O are defined as 

n D = t((T^-a^)-q r xQ p /2V3. (B12) 
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